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SUMMARY 


This paper deals with the coupling of NASTRAN to another finite element 
program which has been developed by the General Motors Research Laboratories 
for the static analysis of automotive structures. The two programs were 
coupled together to use the substructuring capability of the in-house program 
and the normal mode analysis capability of NASTRAN. Modifications had 
to be made to the NASTRAN program in order to make the coupling feasible. 

This information can be of use to other NASTRAN users since there are many 
in-house finite element programs that are specially designed for particular 
problems or have capabilities not found in NASTRAN. By coupling the two 
programs together, the capabilities of both programs can be utilized. 


INTRODUCTION 


An interface program was written to allow an in-house finite element 
program to be used for the static analysis and the NASTRAN program for the 
normal mode analysis of automotive structures. By using the in-house finite 
element program, the following benefits were gained: flexibility of substruc- 

turing, an extended element library, an easily modified program for particular 
problems, and reduced cost of execution. A majority of the geometric data 
had already been prepared and checked in a format compatible with the in- 
house program, thus, the coupling of the two programs saved redefining and 
•debugging the data for the NASTRAN format. Since the in-house finite element 
program does not have dynamic capabilities and the NASTRAN dynamic software 
is reputed to be one of the best, it becomes a logical choice to couple these 
two programs rather than develop a dynamic analysis code or use another 
program. 


This paper describes how the two programs are coupled together. The 
information contained in this paper is based on the use of NASTRAN version 
12.1 on an IBM 370/165 computer. 
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EXPLANATION OF INTERFACE PROGRAM DEVELOPED 


When a substructure is modeled with the in-house code, the grid point 
data and the reduced stiffness matrix for the kept grid points are stored 
in a partitioned data set. When a normal mode analysis of one of these sub- 
structures is desired, this data, along with other user supplied data, is 
loaded into NAS TRAN for the dynamic portion of the analysis by means of the 
interface program. 

The input deck to the interface program is the NASTRAN executive control 
deck, the case control deck, and the bulk data deck, excluding grid point and 
element data. The interface program reads the input deck, converts any 
EBCDIC characters to BCD, and stores the card images in an output file 
for processing by NASTRAN. When processing the bulk data deck, the stiffness 
matrix for the substructure is inserted in its sorted position by placing 
the stiffness coefficients on DMIG cards. Likewise, GRID cards are also 
inserted in the bulk data deck in their correct sorted position. 

An additional program processes the original data used to generate the 
substructure to determine the mass at each grid point by means of the lumped 
mass method. The mass at each grid point is then partitioned out to the 
kept grid points surrounding the grid point. This program then punches C0NM2 
cards for the masses at each of the kept grid points. These cards are then 
read in as part of the bulk data input to the interface program. 

Dummy rod elements are used to define the shape of the structure for 
plotting the mode shapes. To avoid changing the stiffness matrix, the rods 
are given an elastic modulus value of zero. 

Grid point constraints imposed on the structure in the in-house finite 
element program are duplicated in the NASTRAN bulk data deck. This avoided 
calculating extraneous eigenvalues. 

The stiffness values on the DMIG cards are read into the normal mode 
analysis format by means of the ALTER program shown below. 

ALTER 26,27 

MTRXIN, ,MATP00L,EQEXIN, SIL, /STIF, , /V,N,LUSET/V,N,N0M1/C,N,O/C,N,O $ 

SMA1 CSTM ,MPT , ECPT , GP CT , DIT /KGGY , , GPST/ V , N , N0GENL / V ,N ,N0K4GG $ 

ADD KGGY , STIF/KGGX/C ,N , (1 .0 ,0 . 0) /C ,N , (1 .0,0 . 0) $ 

CHKPNT KGGX , GPST $ 

END ALTER 

This is similar to the alter program given in reference 1. 
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MODIFICATIONS TO NASTRAN 


Using the DMIG bulk data cards for inputting the stiffness matrix into 
NASTRAN creates several problems. Even though the stiffness matrix is 
symmetric , the entire matrix has to be read in since the direct matrix input 
method does not take advantage of symmetry. This requires reading in nearly 
twice as many cards as would normally be needed. When reading in large 
matrices, such as a 184 by 184, these extra values consume several minutes of 
computer time in order to pass through the input file processor. 

Another disadvantage of the direct matrix input method is that only 
single precision values for the stiffness matrix can be read in. For the 
IBM 370, single precision is only six significant figures, therefore, round- 
off errors may significantly affect the accuracy of the solution. In order 
to read in six significant figures in an E field format, the large field 
format has to be used. A sorting problem developed when the large field 
format was used on the second DMIG type header card. The interface program 
outputs the DMIG cards in the correct sorted order such that GJ and CJ were 
in an increasing numerical order. However, NASTRAN sorted the deck with 
CJ in a decreasing numerical order. When the cards were input with CJ 
in a decrasing order, NASTRAN sorted them into the correct order of increasing 
CJ number. This problem could not be resolved, so it was reported to the 
NASTRAN system office for further study. 

The sorting problem does not occur if the small field format is used 
for the second type header cards. However, if the small field format is 
used, a stiffness value must be placed on this header card, and the small 
field width does not allow inputting a stiffness value of six significant 
figures in an E field format. To overcome this restriction, a completely 
constrained fictitious grid point with a zero stiffness value is placed in 
this field. 

After the sorting problem was corrected, problems still occurred in 
trying to pass a large number of cards (20,000) through the input file 
processor. Nine minutes of central processing time failed to pass these 
cards through the XS0RT subroutine. The problem was traced to a double 
D0 L00P in the XS0RT subroutine in which a check was made for duplicate 
continuation cards. It appeared that in this subroutine each continuation 
card was compared to every other continuation card, resulting in (N)*(N-l)/2 
comparisons. Approximately 19,000 of the 20,000 cards were continuation 
cards resulting in approximately 1.8 x 10® comparisons. At about 5 micro- 
seconds for each comparison on the IBM 370/165, it would have taken about 
900 seconds to complete this cycle. By removing this double D0 L00P without 
any other modifications to the subroutine, the 20,000 cards were processed 
by the XS0RT subroutine in 138 seconds. 

Since most of the bulk data deck was created by the interface program 
which had been thoroughly checked out, and since the data was already sorted, 
it was felt that most of the checks performed by the XS0RT subroutine were not 
needed. A very brief, modified version of the XS0RT subroutine was created 
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especially for use with the interface program (Appendix) . This modified 
subroutine reads in the data either from direct input, a UMF tape, or a 
check point tape, and outputs the data on the NPTP file tape. This subroutine 
can only be used for UMF tape or check pointed tape problems that require no 
corrections to the bulk data deck. The new subroutine further reduced the 
time to process the 20,000 cards from 138 to 79 seconds. 

Initially, an attempt was made to read in the entire mass matrix by 
means of DMIG cards, using the ALTER program given in reference 1. However, 
the mass matrix MGG must be opened before the program reaches the ALTER DMAP 
instruction where the mass values, read in by the ALTER program, are added to 
MGG. If not, MGG will be an "ill defined matrix," and the addition would 
not be possible. The mass matrix must be opened by some means at DMAP 
operation 28, where MGG is formed. The easiest method is to put in a small 
mass at any grid point. For the problems that were considered only the 
diagonal mass terms were needed, and it was found to be easier to read 
in the mass matrix by means of the C0NM2 cards. This method avoids the 
ill defined mass matrix problem. 


OMITTING ROTATIONAL DEGREES OF FREEDOM 


At first, all the degrees of freedom of each grid point were passed to 
NASTRAN for the dynamic analysis. The rotational degrees of freedom were 
then omitted in NASTRAN, since only the translational degrees of freedom 
were needed to adequately define the fundamental modes of the structures. 

The substructuring program was changed so that it performs a Guyan reduction 
of the stiffness matrix, omitting the rotational degrees of freedom, before 
passing the stiffness matrix to the interface program. This change has 
several advantages. Most of the structures considered are made up from several 
substructures, resulting in stiffness matrices which are densely populated. 

The Guyan reduction process, which tends to fill up the stiffness matrix, had 
little effect on the fullness of the stiffness matrix passed, but it did 
significantly reduce the size of the matrix. The structure which initially 
required 20,000 cards to define the stiffness matrix was reduced to 6,000 
cards with very little change in the eigenvalues and eigenvectors. Since 
the number of cards was greatly reduced, the input file processor time 
was reduced from 79 seconds to 27 seconds. 

This method has the advantage that the reduction process is performed 
with the original sixteen significant figures for each stiffness coefficient, 
instead of the rounded off six significant figures passed to NASTRAN when 
the rotational degrees of freedom are left in. Because these rotational 
degrees of freedom are omitted before being passed to NASTRAN, these degrees 
of freedom have to be constrained in the bulk data deck of the NASTRAN run. 

Back substitution to obtain deflections for the rotational degrees of freedom 
is not possible with this method. 



SPACE FRAME VERIFICATION PROBLEM 


Figure 1 shows a diagram of the structure used to verify the interface 
program. The space frame structure shown consists of 22 grid points 
connected together by 32 rod elements. The space frame was modeled entirely 
in NASTRAN, and the resulting eigenvalues and eigenvectors were compared to 
those obtained by modeling the space frame in the sub structuring program and 
passing the stiffness matrix over to NASTRAN by means of the interface pro- 
gram. A NASTRAN- generated mass matrix was used for both runs, being read in 
as C0NM2 data in the interface case. The values for the 66 eigenvalues and 
eigenvectors agreed to four significant figures. This example problem pro- 
vided an assessment of the loss in accuracy due to the restriction of passing 
only six significant figures for the stiffness values. A comparison of the 
static analysis results between the two finite element programs produced 
exact agreement for the displacements. 

The eigenvalues and eigenvectors for the analysis in which the ro- 
tational degrees of freedom were omitted before the data was passed to 
NASTRAN were closer to the values obtained by using NASTRAN directly than those 
values in which the stiffness values for the rotational degrees of freedom 
were passed to NASTRAN for omitting. 

Omitting the rotational degrees of freedom with the substructuring program 
for the space frame increased the number of card images passed from 1,473 
to 3,674. The reason for the increase was that the stiffness matrix for the 
space frame was very sparsely populated. The Guyan reduction process, 
although it reduced the size of the matrix, produced an almost fully populated 
matrix. 


CONCLUDING REMARKS 


The interface program has been used on considerably more complicated 
structures than the space frame shown. The largest structure analyzed was 
a car body composed of 15 substructures which together contained over 3,000 
degrees of freedom. The substructuring program was used to reduce the struc- 
ture to 186 translational degrees of freedom for passage to NASTRAN for the 
normal mode analysis. The dynamic analysis of the structure was successfully 
completed in 425 seconds with the input file processor requiring 137 seconds. 
The modifications made to the XS0RT subroutine, and the use of the substructure 
program made this analysis technically feasible and computationally economical. 


REFERENCE 


1. Jones, T. C. and Pinson, L. D. , "Adaptation of NASTRAN to the Analysis 

of the Viking Space Vehicle," NASTRAN User's Experience , Sept. 12-15, 
1971, p. 265, Vol. 1. 


115 



o n o n 


APPENDIX 


C^S1RT S READS T RU1 K X CATA CARDS FROM T HP_ INPUT TAPE ADJUST 5 HE 

i: POSITION, AND __ __ 

^tc's THE IfSULT NO SORTED IMAGES ON THE NEW PROBLEM TAPE. 

C * PL&CES B I ^ SKI * 3 IMSK 2 » 8 |.M$K3 ,BI MSK A, BI MSK5 » B IMSK6 

BKVSK't , P.KMSK2," SHIFTS . .. 

CTAR,PLLS, DOLLAR, STARL, SLASH, SFTM, BLANK 
m , r ur >\ BK,BUF,CCNT,CnRSZ, ECHOS, ECHOU, FCNT , OPTP , 

1 COT TAP , PT ST ,P SHI F T , SFT A , ST ARSW ,ORF , ANDF , — — - 

? TRTAl ,T ST , D V F , P ID — ~ 

nflrwlCN H?m"3*<”iElDS?32<,HE4DNS32<,I8LKD4*2<.COCNm<Z 
ol^lcN RK24<^MK^4<tI8UF1320<tI8UF2?20<t ! BUF3?2< *KPARNT,2< 
Cl MENS ICN I BUF l A^2< » IBUF2A%2<*NS0RT??< 


FIELDS, 

TO RIGHT, 


INTEGER 
I ntfcfr 

INTEGER 

INTEGER 


i si . 0l « ♦<.. * «H0£tfas»2Si2^ 

1 f DUM1 , IUEDIT . — 

COMMON /OUT PUT /OUM?%5 fc< » HE AD1%32< , HE AD2%32< , HEA03T?32< — 

COMMON /XMDMSK /DUMA^S <, I CPFLG 

COMMCN/FSCRT/BUFTK — 

CrSm»™MIFSKl*6<,B.«K2«<,^ 

1 BIMSK6 ,BKMSK1X8<,BKMSK2 ’ * F I » c^bTAnk MKATS' 

2 JCCN1 1 ICCN ? , STAR , PLU S, DOLLAR, STARL, SLASH, S FT M, MASK, BLANK, Al_ IS - _ 

CC MM ON /STAPID/ KR AP*12<»KUMF 

ECU! VALENCE lBK*l<,flK'MSKl*5«, *MK«1<»BIMSK2*2<<. 


?MK B 
f SET A 
%MKE 

DATA HE ACU / 1C* AH 
A- H A T , 4H 
DATA HE ACS/ ! 1*AH 
4H D A, AH 


,8IMSK5*1.<<»*INF 
, SH I FTS72<< » TMKD 
,R IMSK5%2<<» -?MKC 


, BT MSK2$1<<* 
,B I MSK 2%2<<_» 
, B IMSKA%1<< 


.AH 
, AH 
t AH 
A, AH 


N, AH 
E , AH 

O, AH 
E, AH 


DATA 

HEACK/ 

?*4F , AH 

f AH 

AH ? 

, AH. * , AH 3 

f AH • • 


AH 6 

, A H . « ,AH 7 

f AH • « 


A Hi 0 

, AH. »5*AH 

/ 

DATA 

CDCNT /AECARC ,AHCOUN, 

AHT / 


0 , AH 
K , AH 
T , AH 
H, AH 
, AH 
* AH 
.AH 


, AH B U, AH 
E, AH C H, AH 
E C, AH B* AH 

0 , 1 0*AH 

, AH . 

A , AH.. , AH 
8 , AH.. , AH 


L K , AH 


J 3 
J Lf 
7 

, AH 1 


. 9*AH 


L * AH K » 


, AH . . , 


,AH 


• • 


.AH 


• • 


/,NS0RT/AHXS0R,A HT _ /. 


DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 


ECHOU, ECHOS /2*0/ 

#1 > A</ 077T777007777 ,07Z77"00007T7T707700^000T777^0^^^ 

m Kfl , V KB , 1 NF , SET A/ 00000007 77777 ^03777^777777.07777777 

mk C /O 0 0 77 7 7 7 77 7 77 / , MKD /07 7 7 777 007 77 7 ( ’ T 

T ENC1, I END? / AH END D, AHA T A / , IEN0/0/ , ISEQ/O/.. ICCgR K/O/ _ — _ 

ST0BtPlU3 »DOLCAR» STARt /4H000*»4H£000» 4H$000»^H^000/ 


I T APE 1 , IT APE2 , I T APE3, IT APE4, 1 T APE 5/ 301 ,302 . 303.3 04, 305 / 

KIN/O/, UMF/AHtJMF / 

C INITIALIZE XSORT 

IPESTR « -I APPRC - 

IF7.KUMF.LE .0< GO TO 90 — — 

KIM' I - - 

CALI OPE N?$ 50 ,UMF , BUF^l < , 2< d . T . n nv*7i<FR~ 

C FIND PARTICULAR BULK DATA FILE ON UMF AS REQUESTED BY _USER^ 


116 



10 CALL R E AD? 1 30, $60, UMF , p I 0 » 1 , 1 , 1 FLG< 

IF?KUMF-PLC< 30,80, 70 
?0 CALL SK PF I L ?U W F , 1 < 

GC TO 10 

30 WR I TESCUTT AP ,35< K!JMF 

35 FnRMATS62H0***USER FATAL MESSAGF 701 , REQUESTED BULK DATA DECK 18, 
1 24H NOT ON USER MASTER F LLE< 

N0G0#-1 

CALL CLOSESUMF, 1< 

RETURN 

50 Ws R I T FSOUT T A P ,5 5< 

55 FORM ATS 1 HO , 2 3X , 2 8H 202 , UMF COULD NOT BE OPENED< 

GC TO 1800 
60 WRITF?0UTTAP,65< 

f 5 F 0 R M A T ? 1 H 0 , 2 3 X ♦ 2 3 H 203, ILLEGAL EOR ON UMF< 

GC TC 1800 

80 CALL CLCSE?UVF,2< 

c ####*# a## ######.## *##«## ##«####### mm mum* *###*« ###«#« 

90 CALL INITCC 

I FST ECHO. L T» C< GC TO 110 
IF? IECHO.EQ. 1< ECHOUtf 1 
IFSTECHO.TQ. ?< ECHOS# 1 
I F? I ECHO.NF. 3< GO TO LOO 
ECFOU# 1 
ECHOS# 1 

100 IF?ICPFLG.NE.O< ECHOS# 1 
110 CONTINUE 

C START WORKING SCRT EUFFER BELOW GTNC I/O BUFFERS 
II# 5* TBUFSZS 1 
IBUFPG# IIG 62 

IBUFLG# COP S Z? BUF?1 < , ICOMQNK - 21 ’ 

IFSIBUFLG-IBUFBG.lt. 210< CALL MESAGES-8 , 0 , NSORT< 

ITAPE# ITAPE1 
JT APE# I TAPE 2 

CALL OPEN ? $ 1750 , NP TP, B UF?4* IBUFSZ £1 < * 3< 

CALL WRITE SNPTP , I BL KO A, 2 , 1< 

IF SIRESTR.GT.O .OR. KIN.GT. 0< GO TO 1400 
190 READ S!NTAPE,190< IBUF1 

CALL XF AD J 1 SI BUFISI < , L SH TFT, 0< 

IF SIBUFISK.FQ. IEND1 . AND. TRUE 1S2C. EQ. I FN02< GO TO 560 
190 FORMAT ?20A4< 

CALL WRITE SNPTP , IBUF1 ,20, 1< 

GO TO 180 

560 CALL FCF SNPTP< 

CALL CLOSE ?NPTP,1< 

RETURN 

1400 IF SKIN.GT.OK GO TO 1430 

CAU OPEN ?*174C,0PTP,BUF?K,0< 

1410 CALL READ SS1730, $1710, OPTP,I8UF3, 2,1, IFLG< 

IF X I BUF3?1< . FO. IBLKCASK. AND. I BUF3S2< • EQ.IBLKDA?2<< GO TO 1420 
CALL SKPFIL SOPT P , £1 < 

GC TC 1410 

1470 CALL READ % $ 1240, $1 7 10, OP TP , I BU FI , 20 , 1 , 1 FLG< 

CALL WRITE SNPTP, IBUF1 ,20, 1< 

GO TO 1420 

1240 CALL CLOSE SQPTP,1< 

GO TO 560 
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1430 CALL OPEN % $ 50 »UMF , BUF?1<, 2< 

CALL PE AD? $1250, $1710,UMF, IBUF1,20, 1,IFLG< 

CALL WRITE ?NPTP , I BUF1 ,20, 1< 

GO TO 1430 

1250 CALL CLOSE ?UMF,1< 

GO TO 560 

1710 WRITE?0UTTAP,1711< 

1711 FORM AT? 1HG , 23X » 2 7H 211, ILLEGAL EOR ON SCRATCH< 

GO TO 1800 

1730 WRITE? OLTTAP, 173 1< 

1731 F0RMAT?1H0,?3X,24H 213, ILLEGAL EOF ON OPTP< 

GO TO 1800 

1740 WRITF?0UTTAP,1741< 

1741 FORMAT? 1H0,23X,25H 214,0PTP COULD NOT BE OPENED< 
GO TO 1800 

1750 WR I TE?OUTT A P , 1 7 5 1< 

1751 FORMATTING, 23X ,29H 215,NPTP COULD NOT BE OPENEO< 
GO TO 1800 

1800 WR I TE?OUTTAP *1801< PLUS 

1801 F0R M AT?A1,2 3H***SYSTEM FATAL ME SS AGE< 

CALL M ES AG E?-37 , 0 , NSORT< 

NCGO # -1 

RFTUPN 

END 
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